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Abstract 

A scale-dependent effective action for gravity is introduced and an exact 
nonperturbative evolution equation is derived which governs its renormal- 
ization group flow. It is invariant under general coordinate transforma- 
tions and satisfies modified BRS Ward-Identities. The evolution equation 
is solved for a simple truncation of the space of actions. In 2 + e dimensions, 
nonperturbative corrections to the /3-function of Newton's constant are de- 
rived and its dependence on the cosmological constant is investigated. In 4 
dimensions, Einstein gravity is found to be "antiscreening" , i.e., Newton's 
constant increases at large distances. 



1 Introduction 



In many of the traditional approaches to quantum gravity the Einstein-Hilbert 
term has been regarded as a fundamental action which should be quantized along 
the same lines as the familiar renormalizable field theories in flat space, such as 
QED for example [|]]. It was soon realized that this program is not only technically 
rather involved but also leads to severe conceptual difficulties. In particular, the 
nonrenormalizability of the theory hampers a meaningful perturbative analysis. 
While this does not rule out the possibility that the theory exists nonperturba- 
tively, not much is known in this direction. However, it could also be argued that 
gravity, as we know it, should not be quantized at all, because Einstein gravity 
is an effective theory which results from quantizing some yet unknown funda- 
mental theory. If so, the Einstein-Hilbert term is an effective action analogous 
to the Heisenberg-Euler action in QED and it should not be compared to the 
"microscopic" action of electrodynamics. 

It seems not unreasonable to assume that the truth lies somewhere between 
those two extreme points of view, i.e., that Einstein gravity is an effective theory 
which is valid near a certain nonzero momentum scale k. This means that it arises 
from the fundamental theory by a "partial quantization" in which only excitations 
with momenta larger than k are integrated out, while those with momenta smaller 
than k are not included. (The interpretation of the Einstein-Hilbert term as a 
fundamental or an ordinary effective action is recovered in the limits k — ■> oo 
and k — > 0, respectively.) An "effective theory at scale k" , when evaluated at 
tree level, should correctly describe all gravitational phenomena which involve 
a typical momentum scale k acting as a physical infrared cutoff. Only if one is 
interested in processes with momenta k! -C k, loop calculations become necessary; 
they amount to integrating out the missing field modes in the momentum interval 
[k',k\. 

We shall regard the scale-dependent action for gravity, henceforth denoted 
Tfc, as a Wilsonian effective action which is obtained from the fundamental ("mi- 
croscopic" ) action S by a kind of coarse-graining analogous to the iterated block- 
spin transformations which are familiar from lattice systems M . In the continuum, 
will be defined in terms of a modified functional integral over e~ s in which the 
contributions of all field modes with momenta smaller than k are suppressed. In 
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this manner Tk interpolates between S (for k — > oo) and the effective action T (for 
k — > 0). The trajectory in the space of all action functionals can be obtained as 
the solution of a certain functional evolution equation, the exact renormalization 
group equation. Its form is independent of the action S under consideration. The 
latter enters via the initial conditions for the renormalization group trajectory; 
it is specified at some UV cutoff scale A: T\ = S. If S is a truly fundamental 
action, A is sent to infinity at the end. The renormalization group equation can 
also be used to evolve effective actions, known at some point A, towards smaller 
scales k < A. In this case A is a fixed, finite scale. In this framework, the 
(non)renormalizability of a theory is seen as a global property of the renormal- 
ization group flow for A — > oo. The evolution equation by itself is perfectly finite 
and well behaved in either case, because it describes only infinitesimal changes of 
the cutoff. 

In this paper we shall give a precise meaning to the notion of a scale- 
dependent gravitational action r^g^] and we shall derive the associated evo- 
lution equation. We employ a formulation in which the metric is the fundamental 
dynamical variable. Alternative approaches based upon the spin-connection and 
the vielbeins are also possible, but they will not be considered here. By using a 
variant of the background gauge technique we are able to make r/Jg^] invariant 
under general coordinate transformations. This property is very important if one 
wants to find nonperturbative solutions of the evolution equations in terms of sim- 
ple truncations of the space of actions. Our construction of Tk [g^u] parallels the 
definition of the "effective average action" [Q, [| which was widely used recently 
H B B SM3 The remarkable successes of this method in flat space are partly due 
to the fact that it allows for nonperturbative solutions when no small expansion 
parameter is available, and that 1^ has a built-in infrared cutoff. Therefore the 
low-momentum behavior of (almost) massless theories can be investigated even 
in cases where IR divergences render standard perturbation theory inapplicable. 
For the purposes of quantum gravity, both of these features are very welcome, 
of course. In fact, in quantum cosmology one of the most intriguing questions is 
how quantized Einstein gravity behaves at extremely large distances. It has been 
argued [|T^, |TJJ that in presence of a nonzero cosmological constant there should 

1 For related work using similar techniques see refs. (I 



11,0. 
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be very strong renormalization effects in the infrared which might even provide a 
mechanism for a dynamical relaxation of the cosmo logical constant. The method 
which we are going to develop would be ideally suited to study problems of this 
type. Since only long distance physics is involved here, there are good chances 
that this can be done without knowing the microscopic theory of quantum gravity. 
(See ref. p| for a related discussion.) 

The "effective average action" used in this paper should not be confused 
with the closely related "average action" which was introduced earlier Ul5| . The 
former obeys a more convenient evolution equation while the latter has a simple 
interpretation in terms of field averages. Their precise relation is explained in 
ref. |tL6|| . The average action has been used in a gravitational context in refs.JTj 



18 1 , but no exact evolution equation was formulated. The evolution of the effective 



average action in a gravitational background was studied in ref. [B3] in the context 
of Liouville field theory. For a review of the effective average action and its 
applicaton to Yang-Mills theory we refer to lpO|| . 

The remaining sections of this paper are organized as follows. In section 2 
we give the definition of I\ and derive the exact, nonperturbative renormalization 
group equation. In section 3 we establish the modified Ward identities satisfied 
by Tfc, and we show that the conventional diffeomorphism Ward identities are 
recovered in the limit k — > 0. In its general form, the evolution equation describes 
a flow on the infinite dimensional space of all action functionals. Approximate 
nonperturbative solutions can be found by truncating the space of actions, i.e., by 
projecting the flow on a finite-dimensional subspace. In section 4 we investigate 
the "Einstein-Hilbert truncation" where only the operators / yfg and / ^fgR are 
retained. In section 5 we determine the resulting scale dependence of Newton's 
constant and of the cosmological constant. As an example, gravity in 2 + e and 
in 4 dimensions is discussed in detail. 



2 The Renormalization Group Equation 

In this section we introduce the effective average action for euclidean quantum 
gravity in d dimensions and we derive the exact renormalization group equation 
which governs its scale dependence. 

We are going to employ the background gauge fixing technique pl[ |22 which 
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means that we decompose the integration variable 7 M ^(x) in the functional integral 
over all metrics according to 

lnv{x) = 9^u(x) + h^(x) (2.1) 

Here g^ is a fixed background metric so that the integration over 7^ may be 
replaced by an integration over h^ v . We consider the following scale-dependent 
modification of the generating functional for the connected Green's functions 



= J Vh^VC^VC, exp {-S[g + h}- S gi [h; g] ( 2 - 2 ) 

-S gh [h, CC; g\ - A k S[h, C, C; g] - S SOUICC } 

Here £[7] = S[g + h] is the classical action which is assumed to be invariant under 
the general coordinate transformations 

= C-vliiv = v p d p >y^ + d^vf-fpv + dyV^^p (2.3) 

where C v denotes the Lie derivative with respect to the vector field v M . For the 
time being let us also assume that S is positive definite. 

Furthermore, S g f denotes the gauge fixing term for the gauge condition 
Fp(g,h) = 0, 

S gi [h; 9} = 7^ J d d x^Tgr u FpF u (2.4) 
and Sgh is the action for the corresponding Faddeev-Popov ghosts C M and C M : 

d d x Cp r v £c {g a p + M (2.5) 

The Faddeev-Popov action is obtained along the same lines as in Yang-Mills 
theory: one applies a gauge transformation 

$hp u £-v1fAU £v(9l!V ~l" hpy) 

(2.6) 

to Fp and replaces the parameters v M by the ghost field C M . The integral over C M 
and Cp provides a representation of the Faddeev-Popov determinant det[5Fp/5v u ] 
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then. In eq. fl2.5|) we introduced the constant (proportional to the Planck mass) 

k = (32ttG)~ 1/2 (2.7) 

where G denotes the bare Newtonian constant. In principle our construction 
works for an arbitrary background gauge fixing. It is particularly convenient to 
use a F p which is linear in the quantum field h pu : 

F„ = V2KFf[g]h aP (2.8) 

We shall mostly employ the harmonic coordinate condition for which T p is the 
following first order differential operator constructed from g pv : 

Tf = 6^3, - l -g aP D, (2.9) 

The covariant derivative D p involves the Christoffel symbols T p of the background 
metric g pv . For the gauge fixing Q2.8p with ( |2"13| ) the ghost action reads 

S gh [h,C,C-g] = ~V2j d*xV§C,M\g,g]» v C u (2.10) 

with the Faddeev-Popov operator 

M[g,gTu = g w 9 uX ^g pv D a + g au D p ) - g?°g» X D x g av D p (2.11) 

The essential piece in eq. ( p.2|) is the IR cutoff for the gravitational field h pv 
and for the ghosts: 

A k S[h,C,C;g] = \k 2 ( d d xjgh pv R%™W vpa h pa 

2 J (2.12) 

+V2jd d x^gC p Rt[g}C» 

The cutoff operators -Rf rav and _Rf h serve the purpose of discriminating between 
high-momentum and low-momentum modes. Eigenmodes of —D 2 with eigenval- 
ues p 2 ^> k 2 are integrated out in ( |2.2| ) without any suppression whereas modes 
with small eigenvalues p 2 <C k 2 are suppressed by a kind of momentum depen- 
dent mass term. The operators -R| rav and i?| h describe the transition from the 
high-momentum to the low-momentum regime. Either of them has the structure 

R k [g] = Z k k 2 R^\-D 2 /k 2 ) (2.13) 
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where the dimensionless function interpolates smoothly between i?(°'(0) = 1 
and lim W°'(u) = 0. A convenient choice is for example 

R(°\u) = u[exp(u) - I}' 1 (2.14) 

The factors 2% are different for the graviton and the ghost cutoff. For the ghost 
Z]~ = Zf 1 is a pure number, whereas for the metric fluctuation Zk = Zf av is 
a tensor constructed from the background metric g^ u . In the simplest case one 
would take 

(Z& w y vpa = g^g ua Zf™ (2.15) 

In section 4 we shall employ a slightly more refined choice. There we shall also 
explain how the factors Zf h and -Zf rav should be choosen. Note that the cutoff 
action (|2.12| ) is quadratic in the quantum fields h^, C 1 and C^. This is an 



important prerequisite for obtaining a tractable evolution equation later on. The 
requirement of a quadratic A^S forces us to use the covariant Laplacian D 2 = 
g^D^Dy in the background metric as the operator which discriminates between 
high-momentum and low-momentum modes. 

In ( p.2|) we coupled h^, and to the sources t^, and cx M , respectively: 

Source = - / d^Mt^h^ + a^ + a^C^ 

J L (2.16) 

+^ U C C (g^ + V) + TpCdvC"} 

The sources (3^ and r M couple to the BRS variations of h^ v and C^ 1 , respectively. 
In fact, it is not difficult verify that S + S g f + S g h is invariant under the BRS 
transformations (e is an anticommuting parameter) 

8 e g»v =0 

(2.17) 

8 e C» = eK~ 2 C p d v C» 
5 £ Cf, = e (an)' 1 Ff, 

Given the functional W^, we introduce fc-dependent classical fields 



6 



and we formally solve for the sources (t^ , cr M , a^) as functionals of the fields 
(fyitK > ^ i £n) an d °f (P' 11 ' j T fi j Then the Legendre transform l\ of W*. de- 

pends on the classical fields and parametrically on (3, r and g: 

r k [h,Z,Z;l3,T ] g] = j d d x^{n, v + a^ + a%}-W k %a 1 a-f3 1 r-g\ (2.19) 

By definition, the effective average action T k obtains from T k by subtracting the 
cutoff action A k S with the classical fields inserted: 

r fc [7i, £, £; A r; = f fc [/i, £, 6 ft r; g] - A k S[h, £, £; (2.20) 

It is convenient to define the metric 

9fiu{x) = g^(x) + hp,(x) (2.21) 

as the classical analogue of the quantum metric 7^ = g^ + and to consider 
functional of g^ rather than h^ u : 

^k\9iu>,9^,^,^P,r] =^k[g^-9^^' x ,^;P,r;g IJ , v ) (2.22) 

The main virtue of the background technique employed here is that the functional 
Tfc is invariant under general coordinate transformations where all its arguments 
transform as tensors of the corresponding rank: 

r fc [$ + £,$] = r fc [$] , ^ = {g,u,g^,e,^f3^,r,} (2.23) 



Note that in Q2.23 ), contrary to the "gauge transformation" (2.6), also the back- 



ground metric transforms as an ordinary tensor field: Sg^ = C v g^ u . Eq. ( [2.23 
a consequence of 



is 



W k [J + C v J} = W k [J] , J = {t^, <7 M , cx^; r M ; g^ u } (2.24) 



This invarianc property follows from (|2.2f ) if one performs a compensating trans- 
formation on the integration variables h^, C 1 and C^. At this point we assume 
that the measure is diffeomorphism invariant. 

The general coordinate invariance of T k is of major practical importance be- 
cause if we know a priori that no symmetry-violating terms are generated during 
the evolution it is sufficient to use truncations which consist of invariant combina- 
tions of the fields only. The conventionally defined effective action of the metric, 
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r[fiv], obtains in the limit of a vanishing IR cutoff by setting the ghosts, (3 and r 
to zero and by identifying g^ v with g^ u : 



r bv] = Jim T k [g^ v , g^,0,0] 0,0] 



(2.25) 



As a consequence, r[g^] is invariant under 5g^ v = C v g^ v . Even though we are 
mostly interested in the functional 



Tk[g^v] = F k [g^ u , g^, 0, 0; 0, 0] 



(2.26) 



which depends on g^ v only, an exact renormalization group equation can be for- 
mulated only if one keeps track of the dependence on £, £ and g as well. For the 
derivation of the (modified) BRS Ward identities satisfied by Y k the dependence 
on p and r must be retained in addition. 

The derivation of the evolution equation for T k proceeds as follows. Taking 



a derivative of the functional integral ( |2.2| ) with respect to the renormalization 
group "time" t = In k one obtains, in matrix notation 



- d t W k = -Tr 



<h®h> (d t R k y 



)hh. 


- Tr 







< C ® C > (d t Rk) 



(2.27) 



Here R k is a matrix in field space whose non-zero entries are 



( 



Ri 



\ fiupa 
'Jhh 



k 1 {RT v W vpa 



(2.28) 



(Jfc) a = V2Rf[g] 



The RHS of ( 2.27| ) can be expressed in terms of T k by noting that the connected 
two-point function 



Gij(x,y) = < Xi{x)Xj{y) > -<Pi(x)(pj(y) 
1 S 2 W k 



(2.29) 
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x)g{y) SJ i (x)5J^(y) 



and 



x,y) 



(2.30) 
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are inverse matrices in the sense that 

id„, r< t„ ».\ -p( 2 )i^ 



J d d y^)G ij (x,y)f ( £ )jl (y,z) = 8 l i 



i S{x-z) 



(2.31) 



0t* 



Here we used the shorthand notation \% = = an d y?j = 

{/i, £, £}. Thus one obtains the evolution equation 



d t r k [h,£,t,f3,T;g} = -Tr 



-Tr 



( r * )+ ^)^(^- 



hh 



(2.32) 



If one evaluates the RHS of this equation in terms of position-space matrix el- 

(2) 

ements then 1^ is defined by a formula similar to ( |2.3(J| ) and the integration 
implied by "Tr" has to be interpreted as / d d x \Jg(x). The matrix elements in the 
ghost sector are defined in terms of left derivatives, e.g. 

1 5 1 5T k 



rl 2) V 



k J tth* ~ ^{xj s^(x) 5Uv) 



(2.33) 



For any cutoff which is qualitatively similar to (p.!4j) the traces on the RHS 
of eq. (|2.32| ) are well convergent, both in the IR and the UV. By virtue of the 
factor d t R k , the dominant contributions come from a narrow band of generalized 
momenta centered around k. Large momenta are exponentially suppressed. 

Solving the evolution equation ( |2.32| ) with the appropriate initial condition 
at the UV cutoff scale A — > oo is tantamount to computing the original functional 
integral (|2.2|) . In order to determine the correct initial value I\ we consider the 
following integral equation satisfied by IV 



exp{-r k [h,€,bP,T;g]} 



I 



VhVCVC exp 



S[h,C,C;P,r;g] 



(2.34) 



• exp {-A k S[h -h,C-£,C-&. 
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Here 



S = S + Spf + s. 



d d x VI {P^Ccig^ + V) + ^CUC} (2.35) 



is expressed in terms of the "microscopic" fields (h,C,C). Eq. ( |2.34| ) obtains by 



inserting the definition of Tk into (|2.2) and using 



5Tj 



-VI^ 



-VI ^ (2-36) 



The crucial observation is that for k —>■ oo the last exponential in ( |2.34|) becomes 
proportional to a 5-functional which equates the quantum fields (h, C, C) to their 
classical counterparts: 



(2.37) 



As a consequence, the effective average action at the UV cutoff readsQ 



FA[h,£,£;f3,T;g] = S[g + h] + S g{ [h; g] + S gh [h, f, f; g] 

(2.38) 



-fd d x^g- {r^{g»» + Kv) + T,ed„e} 



It is this action I\ which has to be used as the initial condition for the evolution 
equation. We note that at the level of the functional Th[g] eq.( [2.38j ) boils down to 



f aW = S[ 9fiU ) (2.39) 

As T k involves derivatives with respect to g^ u at fixed g^ v it is clear that the 
evolution equation cannot be formulated in terms of T k alone, however. 

Up to now we assumed that the fundamental action S is positive definite and 
the euclidean functional integral ( |2.2| ) makes sense as it stands. It is well known 
that this is not the case for the Einstein-Hilbert action, for example, because the 
conformal factor has a "wrong sign" kinetic term. Clearly it would be desirable 
to have an evolution equation which can be applied in such CclSGS clS well. It is 



quite remarkable therefore that the renormalization group equation ( 2.32|) , with a 



2 Strictly speaking (2.38) is correct only up to local terms which at most change the bare 
parameters in S. Because the value of the bare parameters has anyhow no physical significance 
we ignore these terms here. 
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properly chosen cutoff, is well-defined even if S and Tk are not positive definite. 
To see this, let us look at the first trace on the RHS of ( |2.32|) and let us concentrate 
on the contribution of a fixed mode <ft contained in the metric. We assume that <p 
is an eigenfunction of with eigenvalue z k p 2 where p 2 is a positive eigenvalue 
of some covariant kinetic operator, typically of the form — D 2 + R-terms. For 
theories with S > 0, the wave function renormalization z k is positive (at least 
for large k). In this case the general rule || || is to define the constant Z k in 
the cutoff Rk, eq. ( [2.131 ), as ^fc — z h because this guarantees that for the low- 
momentum modes the effective inverse propagators r^ 2 ^ + Rk becomes Zk(p 2 + k 2 ), 
as it should be. 

The important question is how Z k should be chosen if Zk is negative. If we 
continue to use Z k = z k , the evolution equation is perfectly well defined because 
the inverse propagator — \z k \(p 2 + k 2 ) never vanishes, and the traces of ( |2.32| ) are 
not suffering from any IR problems. In fact, if we write down the perturbative 
expansion for the functional trace, for instance, it is clear that all propagators are 
correctly cut off in the IR, and that loop momenta smaller than k are suppressed. 
On the other hand, if we set Z k = —z k , then — \z k \{p 2 — k 2 ) introduces a spurios 
singularity at p 2 = k 2 , and the cutoff fails to make the theory IR finite in this 
case. 

At first sight the choice Z k = —z k might have appeared more natural because 
only if Z k > the factor exp (— A k S) ~ exp (— / R k (p 2 ) is a damped exponential 
which suppresses the low momentum modes in the usual way. In this paper we 
shall nevertheless adopt the rule Z = z k for either sign of z k . We shall see that 
at least for the Einstein-Hilbert truncation of section 4 the evolution equations 
are well defined and consistent even though it is difficult to give a meaning to 
the functional integral itself. In the case Z k = z k < the factor exp (+ / \R k \<p 2 ) 
unavoidably becomes a growing exponential and it might seem that this enhances 
rather than suppresses the low momentum modes. However, as suggested by the 
perturbative argument above, this conclusion is too naive probably. Moreover, if 
one invokes the usual prescription of rotating the contour of integration over so 
that it is parallel to the imaginary axis, both the kinetic term and the cutoff lead 
to damped exponentials. 

Furthermore, it is important to note that the constructions in this section can 
be repeated for metrics on Lorentzian spacetimes. Then one deals with oscillating 
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exponentials e lS , and for arguments like the one leading to eq. (|2.37|) one has 



to employ the Riemann-Lebesgue lemma. Apart from the obvious substitutions 
T k — ► — iT k , Rk —iRk, the evolution equation remains unaltered. For Z k = z k 
it has all the desired features, and z k < seems not to pose any special problem. 

3 Modified Ward Identities and 
Consistent Truncations 

We mentioned already that the classical action plus the gauge fixing and ghost 
terms are invariant under the BRS transformations ( |2.17| ). Therefore the BRS 



variation of the total action S to t = S + S g f + S g h + AkS + S sonrces receives contribu- 
tions only from the cutoff and the source terms. If we apply a BRS transformation 
to the integral defining Wk and assume that the measure is invariant we obtain 

< ^sources + 5 £ A k S >= (3.1) 

where 

< O >= e~ Wk J VhVCVC O e~ Stot (3.2) 



Our goal is to convert (3JJ to a statement about the average action IV Because 



the BRS transformation ( 2.17Q is off-shell nilpotent when acting on h^ v and on 



(but not on C M ) one has 

^sources = -EK~ 2 / A Vjf^^tV + V) 

J (3.3) 

If we take the expectation value of ( |3.3| ) and express Wk in terms of T k we find 

r q . e f d 1 j ST' k 5T' k 5T' k 5T' k \ e~ 
< S e S Bomce >= _ / dx -j— I _ _ + - - J + -Y k (3.4) 



K 2 



with 

1 (8A k S 5T' h 5A k S 5T[ 



Here we defined 

T' k = T k -S gi [h;g} (3.6) 
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and we exploited the equation of motion < SStot/SC/j, >= which can be cast in 
the form 



1 



6 



8 



r k [h,£,&P,T]g] = (3.7) 



The variation of the cutoff terms gives rise to 

e 



< 5 e A k S > 



(3.8) 



with 



K = /t 2 Tr 



-\/2Tr 



(3.9) 



-2crVTr 



where (p = {/i, is summed over. From 
identities in their final form: 



and (13.81) we obtain the Ward 



d d x- 



ST' k ST' k 5T' k 5T' k 



(3.10) 



Eq.( p.lOj ) has to be compared to the ordinary gravitational Ward identities [ffijfl 
which are similar to (|3.10|) but with a vanishing RHS. In fact, the contribution Y k is 
due to the cutoff and therefore it vanishes in the limit k — > because R k ~ k 2 — > 
in this limit. Hence the standard effective action IfinF/;, is guaranteed to obey its 
usual Ward identities, and BRS invariance is restored for k — > 0. 

Because the Ward identity ( |3.10|) is derived from the same functional integral 
as the evolution equation, it is automatically satisfied for the exact solution of the 
evolution equation. For approximate solutions of the evolution equation their 
consistency with the Ward identity is not guaranteed, and one may even use 
(|3.10| ) to judge the quality of the approximation ^ , [Lj| . 

The most important strategy for finding approximate (but still nonpertur- 
bative) solutions to the evolution equation is to truncate the space of action 
functionals. Typically one works on a finite-dimensional subspace parametrized 
by only a few generalized couplings. As a first step towards such a truncation one 
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can try to neglect the evolution of the ghost action. This amounts to making an 
ansatz of the following form: 



rfc[#,#,£,£;/3,T] = T k [g] +t k [g,g] + S gi [g-g;g] + S gh {g - g,£,£; g] 

(3.11) 

-Jd d x^g {/TT^ + T^dyt*} 
In ( |3.11| ) we pulled out the classical S g { and 5* g h from T k , and also the coupling 



to the BRS variations has the same form as in the bare action. The remaining 
functional depends on both g^ u and g^ u . It is further decomposed as f k + Y k 
where T k is defined as in ( |2.26| ) and T k contains the deviations for g ^ g. Hence 
by definition 

f k [g l9 } = (3.12) 

T k can be viewed as a quantum correction the gauge fixing term which also van- 
ishes for g = g. The ansatz (|3.11|) satisfies the initial condition ( 2.38Q if 



Ta = S , f A = (3.13) 

and it satisfies the quantum equation of motion ( |3T7| ) exactly. Eq. ( |3.13| ) suggests 
to set T k = for all k in a first approximation. In this case it can be checked 
that if the ansatz (|3.11|) is inserted into the Ward identity ( |3.10|) its LHS vanishes 



identically. Including T k the Ward identity assumes the form 

[d d xC^ 6 -^4 = -Y k (3.14) 

We see that Y k = is a good approximation provided we may neglect Y k . The 
traces which define Y k amount to loop integrals, and if we think in terms of a 
loop expansion Y k is certainly a higher loop effect and may be neglected in a 
first approximation. At the nonperturbative level one can still try to set 1^ = 
and investigate the consequences in concrete examples. In Yang-Mills theory the 
analogous truncation has led to rather encouraging results already ||, [| ||. In 
the next section we shall perform an explicit calculation in this approximation. 

If one inserts the ansatz ( |3.11| ) into the evolution equation ( |2.32| ) one finds 
the following equation for the evolution of T k in the subspace spanned by the 
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ansatz: 



d t Fk[g,g] = 2 Tr 



-Tr 



(-M[9,g) + Rt^y 1 d t Rt[ 



This equation is written down in terms of 



r k [g,g) = T k \g,g,0,0;0,0] 



T k [g} + S gi [g - g;g] +f k [g,g] 



(3.15) 



(3.16) 



r^. is the Hessian of r^^, g\ with respect to g^ v at fixed g^ v . For the harmonic 
coordinate condition, the classical kinetic term of the ghosts, Ai, is given by 
eq.fETLTD. 



4 The Einstein— Hilbert Truncation 



In this section we illustrate the use of eq.( ij.l5| ) by means of a simple example. At 



the UV scale A we start from the classical Einstein-Hilbert action in d dimensions, 

s = itc I ddx ^{~ R ^ + 2A I ' 

and we evolve it down to smaller scales k < A. For the time being we shall not 
try to send A to infinity, so the nonrenormalizability of the theory is not an issue 
here. We are going to use a truncation which replaces in (|4.1| ) the bare Newton 
constant G and the bare cosmological constant A by fc-dependent functions 

G k = Z^k 1 G (4.2) 

and Afc, respectively: 



r k [g,g] = 2n 2 Z Nk J d d x y/g {-R{g) + 2X k } 

Wz Nk J Avff (rfgap) (Ffgp*) 



(4.3) 
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This ansatz is of the form Q3.16 ) with T k neglected and the classical gauge fixing 
term given by with and a = l/Z Nk . (Note that Tf g aj3 = T a fh a p 
because D^g^ = 0.) In order to determine the functions Z^ k and X k we have 
to project the evolution equation on the space spanned by the operators ^fg and 
,JgR. After having inserted the ansatz into the evolution equation we may set 
g^v — g^iu so that the gauge fixing term in (|4.3| ) vanishes. The LHS of the evolution 
equation reads then 



d t T k [g, g] = 2k 2 f d d x y/g \-R{g)d t Z Nk + 2d t (z Nk X 



(4.4) 



On the RHS of ( |3.15|) we have to perform a derivative expansion and retain only 
the terms proportional to J y/g and / yfgR. Equating the result to ( [4.4[ ) we can 
read off the system of ordinary differential equations for Z^ k and X k . They have 
to be solved subject to the initial conditions Z^k = 1 and Aa — A. In this manner 
the renormalization group flow in the space of all action functionals is projected 
onto the 2-dimensional subspace parametrized by G and A. 

In the evolution equation we need the second functional derivative of T k [g, g] 
at fixed g^y. We expand 



r k [g + h,g] = r k [g,g] + 0(h) + rr>; g] + 0(h 

and we find for the piece which is quadratic in h^: 



rr d [/>; g] = z Nk K 2 I d d x jg V -K»»D 2 + u»% 



pa 



(4.5) 



(4.6) 



Here indices are raised and lowered with g pu , and the tensors K and U are given 
by 



1 r 



pa 



+ 6w - r u g P . 



(4.7) 



and 



pa 



5$% + 5*% - r u 9po\ (R - 2A fc ) + - [g^Rpa + g P aR» u 



5»R\ + b»R% + + 5 v a R» p 



Ry a + RY P 



In eq. (|4.8| ) all geometrical quantities are constructed from the background metric.^ 
In order to partially diagonalize the quadratic form fl4.6| ) we write as the sum 



We use the conventions R", 



ppv 



-d v Y 



pp 



Rp,v = R a p<jv an d R = Q^Rpv 
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of a traceless tensor h^ v and a trace part involving <ft = g^h^: 

hf, u = h^ u + d^g^cj) , g^h^ = 
As a consequence, eq.(|4.6|) becomes 



(4.9) 



rr d [kg] 



Z Nk K l I d d x - V [-D - 2A fc + R 



Ad 



-D - 2A* 



d-A 



d 



-R 



(4.10) 



d-A 



d 

The equations for Z^ k and obtain by comparing the coefficients of / y/g and 
/ y/gR on both sides of the evolution equation at g^ v = g^ u . For this purpose we 
may insert an arbitrary family of metrics g^ v which is general enough to identify 
the terms / yfg and / y/gR and to distinguish them from higher order terms in 
the derivative expansion, such as / ^fgR? or / ^/gR^D^D^R, for instance. We 
exploit this freedom by assuming that g^ u corresponds to a maximally symmetric 
space, i.e., that 

j^y [d^pQua g^o-gup] R 



R 



pvpcj 



(4.11) 



R 



pu 



d 



gpvR 



From now on the curvature scalar R parametrizes the family of metrics inserted, 
and it should be regarded as an externally prescribed number rather than a func- 
tional of the metric. For a maximally symmetric background the quadratic action 
boils down to 



rr d [M 



Z Nk K 2 I d d x{h 



pv 



-D - 2\ k + C T R 



'd-2" 



2d 



cf) -D 2 - 2X k + C S R 



with 



Ct = 



d(d- 3) +4 
d(d- 1) 



C, 



d-A 
d 



(4.12) 



(4.13) 
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Before continuing we have to specify the precise form of the cutoff operators 
i?| rav and i?| h to be used in the evolution equation (|3.15|) . Both of them have 



the structure ( [2.13| ) whereby Z k should be adjusted in such a way that for every 
low-momentum mode the cutoff combines with the kinetic term of this mode to 
— D 2 + k 2 times a constant. Looking at ( |4.12| ) we see that the respective kinetic 
terms for and differ by a factor of — (d — 2)/ 2d. This suggests the following 
choice: 



{I -Pi 



\fiupa 



d-2 
~2A 



'-P, 



/ivpa 



(4.14) 



Here 



(P<f),r = (4.15) 
is the projector on the trace part of the metric. For the traceless tensor ( |4.14j ) 



coincides with ( 2.15 ) for Zf rav = Zj^k, and for <p the different relative normalization 
is taken into account. Thus we obtain in the h and the 0-sector, respectively: 



^- 2 T^[g,g] + R 
(n- 2 T^[g,g] + R 



ZfJh 



-D 2 + k 2 R i0 \-D 2 /k 2 ) - 2\ k + C T R 



grav 
k 



d-2 



2d 



, sk —D 2 + k 2 R^\-D 2 /k 2 ) - 2X k + C S R 

(4.16) 

From now on we may set g = g and we omit the bars from the metric and the 
curvature. 

The last missing ingredient for the evolution equation is the Faddeev-Popov 
operator. From ( |2.11 ) one obtains at g = g 



with 



-D 2 + C V R 



(4.17) 



(4.18) 



In the second part of ( 4.17|) we used (|4.11 ) for a maximally symmetric background. 
Since we did not take into account any renormalization effects in the ghost action 



we set Zf 1 = 1 in i?| h and obtain 



M + Rf = —D 2 + k 2 R^\-D 2 /k 2 ) + C V R 



(4.19) 
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Let us write Sk(R) for the RHS of the renormalization group equation (|3.15 ) 
with g = g. Inserting (|4.16| ) and (|4.19 ) there we arrive at 



with 



S k (R) = Tr T W(A + CtR)' 1 } + Tr 5 \N{A + C s R)~ l 



2Tr y Wo(Ao + C v RY l 



A = _ jD 2 + A ,2 jR (0) ( _ jD 2 /A; 2 ) _ 2 Y fc 

M = (2Z Nk )- 1 d t \z Nk k 2 R^(-D 2 /k 2 ) 



1 - 7flN{k) 



(4.20) 



(4.21) 



k 2 R^{-D 2 /k 2 ) + D 2 R^'{-D 2 /k 2 ) 
where a prime denotes the derivative with respect to the argument and 



rj N (k) = -d t lnZ 



Nk 



(4.22) 



is the anomalous dimension of the operator ^fg~R. The operators and Aq are 
defined similarly to ( ggg ) but with A = and Z Nk = 1, i.e., rj N {k) = 0. Eq.(^2"D|) 
involves traces of functions of the covariant Laplacian D 2 = g^ u D^D v acting on 
traceless symmetric tensors ("T"), scalars ("S"') and vectors ("V™). Because we 
need only the zeroth and the first order in the curvature scalar we can expand 



S k {R) = Tit 



MA 



-i 



+ Trc 



NA 



-1 



2Tr 



v 



- R (C T Ti T [UA- 2 \ + C s Tt s [MA 



NoAq 1 



2C v Tt v 



(4.23) 



M Af 



+0(R 2 



The traces in ( |4.23 ) can be evaluated by taking advantage of the heat kernel 
expansion 

d/2 



Tr 



-±-J~tr{I) j d d x^\l- l -isR + 0{R 2 ) 



(4.24) 



Here I denotes the unit matrix of the space of fields on which D 2 acts. Hence 
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tr(J) is the number of independent field components and in particular 



tr s (7) = f 
try (7) = d 

tr r (J) = i(d-l)(d + 2) 



(4.25) 



Considering an arbitrary function W with a Fourier transform W, the expansion 
of the trace 



Tr{W(—D 2 )} = / dsW(s)Tr 



-isD 2 



(4.26) 



is given by 



Tr[W(-D 2 )} = (4n)- d /hr(I)^Q d/2 [W} J d d x y/g 



+ \Qd,2-i[W} J d d x^R + 0(R 2 



(4.27) 



with 

/oo 
ds(-is) n W(s) (4.28) 

Reexpressing ( |4.2(j| ) in terms of W leads to the Mellin transform (n > 0) 



Qo[W] = W(0) 



(4.29) 



Qn[W] 



Tin) Jo 



dzz n ~ L W(z) 



The next step is to use ( |4.27| ) in order to evaluate ( |4.23|) and to combine S(R) 
with the LHS of the evolution equation, eq. ( t4.4|) . From the coefficients of / ^fg 
we can read off the following equation 



d t (z Nk X k ) = (4 K 2 )- 1 (47r)- d / 2 \tr T (I)Q d/2 [jV/A] 



(4.30) 



+tr s (I)Q d/2 [/V/A} - 2ti v (I)Q d/2 [f\f /A } 
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Likewise / gives rise to 



d t Z Nk = -(12^)- i (4 7 r 



2 \- l (A^r\- d / 2 



tr T (/) {Qd/2-xW/A ~ QC T Q d/2 [Af/A 2 ]} 



+ 



tr s (I) {Q d /2-iW/A] - QC s Q d/2 [M/A 2 ]} ( 4 - 31 ) 



-2tr v (/) {Q d/2 _i[JVo/A)] - 6C v g d/2 [jV /^]} 



In (|4.30| ) and ([4.31 ), A/" and A are considered c-number functions of z which 
replaces — D 2 in (|4.21 ). For every cutoff we define the functions (p — 1, 2, . . .) 



i 



r(n) Jo 

1 /"OO 



r(n) Jo 



n _ lJ R( )(^)- 2/j R(°)'^) 



#°>(z) 

> + i?(°)(2)+w]P 



(4.32) 



for n > 0, ancf] 



$g( w ) = $g( w ) = (l + w )-P 
In terms of the $'s, eq. ( |4.30| ) assumes the form 

d t (z Nk \ k ) = (16/ t 2 )- 1 (4vr)- d / 2 A;42 C /( C /+l)<I>i /2 (-2A fc /A; 2 )-8^ /2 (0) 



(4.33) 



(4.34) 



-d(d+l) VN $ d/2 (-2\ k /k 2 ) 



and (|4.31|) becomes 



d t Z Nk = -(24 K 2 )- 1 (47r)- d / 2 A; 



d-2 



d(d + l) ^ /2 _ 1 (-2A fc //c 2 ) - -^$^^(-2^/^ 



-6d(d - 1) $ 2 /2 (-2A fc /fc 2 ) - -^ /2 (-2A fe /fc 2 ) 



-Ad& d/2 _M -24$ 2 /2 (0) 



Actually eq.(4.33) follows from (4.32) in the limit n\0 



(4.35) 
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Let us introduce the dimensionless, renormalized Newton constant 

9k = k d 2 G k = k d 2 Z Nk 1 G 
and the dimensionless cosmological constant 

X k = k 2 X k 



(4.36) 



(4.37) 



Here G k = Z^ k 1 G is the dimensionful renormalized Newton constant at scale k. 



The evolution of g k is governed by the equation 



d t 9k = [d- 2 + r} N (k)]g k 



From (|4.35 ) we obtain for the anomalous dimension r]N{k): 

VN(k) = g k Bi(X k ) + rj N (k)g k B 2 {\ k ) 



with 



1 



B x (X k ) = -(4tt) 



l-d/2 



d(d + l)^ /2 _ 1 (-2A fc ) - 6d(d - l)$* /2 (-2A fc ) 



d/2\ 



-4d$i /2 _ 1 (0)-24^ /2 (0) 



d/2\ 



B 2 (X h 







d{d+l)& d p_ x {-2\ k )-§d(d-l)Q 



d/2\ 



-2A* 



We can solve ( |4.39| ) for the anomalous dimension in terms of g k and A^: 

gkBi(Xh) 
VN ~ l-g k B 2 (X k ) 

The scale derivative of X k is related to (|4.34j) according to 

d t \k = -(2 - r] N )X k + 32Tcg k K 2 k~ d d t (Z Nk X k ) 

so that 



d t X k = -{2-r lN )X k + ^g k {A'K) l - d l 2 - 



2d(d + l)$^ /2 (-2A fc ) - 8d& d/2 {0) - d(d + l)w$i /2 (-2A 



d/2\ 



(4.38) 



(4.39) 



(4.40) 



(4.41) 



(4.42) 



(4.43) 
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Eqs. (|4.38 ) and ( 4.43|) with Q4.41] ) is the set of differential equations we wanted 



to derive. Once the initial values g& and Aa are given, it determines the value 
of the running Newton's constant and cosmological constant at any scale k < 
A. Although they were derived from a relatively simple truncation, the above 
evolution equations encapsulate nonperturbative effects which go beyond a simple 
one-loop calculation. This is particularly obvious if one expands for instance 
(|4.41|) for small values of g k : 

Vn = g k Bi{\ k ) [l + g k B 2 (\k) + glBl(\ k ) + ■■■] (4.44) 

We observe that 7]^ receives contributions from arbitrarily high orders of pertur- 
bation theory. 

5 Running Newton's Constant and 
Cosmological Constant 

5.1 Near two dimensions 

In d = 2 dimensions / JgR is a topological invariant proportional to the Euler 
number and the quantum theory under consideration has at most finitely many 
(topological) degrees of freedom. In d = 2 + e dimensions, on the other hand, one 



finds a dynamically nontrivial theory with a nonzero /3-function for g k p4, |25 



d t gk = [e + t) N ] g k (5.1) 

Gravity in 2 + e dimensions provides an interesting laboratory for a first test of the 
evolution equation because here the conformal factor of the metric can have both 
a conventional (e < 0) and a "wrong-sign" (e > 0) kinetic term, see eq.( |4.12 ). 
The anomalous dimension has a power series expansion 

Vn = V ( n ] + V$e + 4 ] e 2 + -.. (5.2) 

and therefore 

d t g k = \(l + V $) e + g k + 0(e 2 ) (5.3) 
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Expanding the functions fl4.4U| ) as = b[°1 + B^e + . . . one has 



(o) 5 , fc-Bi° ) 



1 - QkB^ 



(D 9kB? glBTB± 



(5.4) 



The lowest order terms are 



i?i 0) (A fe ) = 2 (l-2A fc )- 1 -4$?(-2A fc )-y 



Bi°\\ k ) = 2$2(-2A fe )-(l-2A fc )- 1 



(5.5) 



We remark that for vanishing cosmological constant, B^ is a universal quantity, 
i.e., it does not depend on the precise form of R^: 

B?\0) = - 3 j (5.6) 

The reason is that the integrand in the integral representation of $f(0) equals the 
derivative of z(z + R^(z))^ 1 ; hence it is sufficient to know that is bounded 
everywhere in order to establish that 

$?(0) = 1 (5.7) 

Unlike $i(0), $f(Aj.) is sensitive to the shape of R^ even for X k = 0. In order 
to be more explicit we evaluate ( |5.5[) at A ^ for the constant cutoff function 
R(°'(z) = 1. Though it does not vanish for z — > oo, it yields at least qualitatively 
correct results || [jj as long as it does not introduce UV divergences into the 
integral under consideration. For $^ and $^ this is not the case and one finds 

$l(w) = $l(w) = (1 + w)- 1 (5.8) 

so that 

B[°\\ k ) = -2(l-2A fc )- 1 -f 

3 (5.9) 

Bi°\\ k ) = (l-2X k y 1 
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As a consequence, we obtain the following answer for the anomalous dimension: 

(o) 38 1 - f§A fc /- in x 

W = ~^r9k- ox (5.10) 



Z a *l-g k -2\ 



k 



Eq.( p.l(J| ) improves on earlier results in refs.^j, [25], 26|. It takes into account 
partially resummed higher loop effects (higher powers of g k ) and it includes the 
effect of the running cosmological constant. 

One of the interesting features of Einstein-Hilbert gravity in 2+e dimensions 
is that the evolution of Newton's constant is governed by a fixed point g* at which 
the /^-function (|5.3j) vanishes. To lowest order in e it is given by 

g * = -eBf\\ k )- 1 (5.11) 

The A-dependence of g* is non-universal. For = 1 we obtain 

3 l-2\ k 

s - = Ts £ T^ (5 - 12) 

Eq.( |5.12j ) is reliable for A,t C 1. In this regime the fixed point g* is UV stable 
if e > and it is IR stable for e < 0. For e > and A^ = this fixed point 
was discussed by Weinberg in the context of the asymptotic safety scenario 



for quantum gravity. Our result for the dependence of g k on the cosmological 
constant can only be obtained in a framework with a proper infrared regularization 
because we are investigating the influence of the relevant dimension-two operator 
on a marginal coupling. (In a sense, the role played by the running cosmological 
constant is similar to the quadratic mass renormalization in four dimensional 
scalar theories.) For e > the theory is asymptotically free. Near the fixed point 
the dimensionful Newton constant G k = g*/k £ vanishes for k — > oo. 

The evolution of A& itself is governed by eq.Q4.43p. For g k ~ g* where g k and 



rjN are of order e, one finds that also the /9-function of A has a zero of order e: 

A * = ~^ M£ (5 - 13) 

This fixed point of the A-evolution is UV stable for either sign of e. We conclude 
that to first order in e and for e > the combined (A, g)-system has an UV stable 
fixed point given by ( |5.13p together with g* = (3/38)e. 
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5.2 Four dimensions 



In d = 4 dimensions, the running of Newton's constant is governed by the following 
functions of the cosmological constant: 

1 

~3% 



B x {\) 



18^(-2A) - 5$J(-2A) + 6$^(0) + 4$}(0) 



B 2 {\) 



1 

6n 



18$|(-2A) - 5$J(-2A) 



(5.14) 
(5.15) 



The dimensionful quantity G k evolves according to 

dtGk = VNGk 



(5.16) 



with the anomalous dimension given by ( [4.411 ). In order to get a feeling for 
the behavior of Gk, let us restrict our attention to the lowest order in g k which 
amounts to keeping only the first nontrivial correction of the expansion in Gk 2 . 
Then 77^ = Si(A fc )^ + . . ., or with g k = k 2 G k = k 2 G + 0(G 2 ), 



r ]N = B 1 (\ k )Gk 2 + 0(G 2 



(5.17) 



First we consider the case where the cosmological constant is much smaller than 
k 2 . Then we may approximate A& ~ in (|5.17|) , and (|5.16[) has the solution 



Gb — Gr 



1-uGk 2 + 0(G 2 k 4 ' 



Here 



1 

6n 



24^(0) -$i(0) 



(5.18) 



(5.19) 



is a pure number, which depends on the function , however. For the exponen- 
tial cutoff ( gTp we have $J(0) = tt 2 /6 and $ 2 (0) = 1, so that 



7T 



CO 



4 / 

7T V 1 I- 1 



> 



(5.20) 



For different cutoff functions the numerical value of u will be slightly different 
but it will still be positive. Therefore eq. (|5.18 ) tells us that Newton's constant 
decreases as k 2 increases; it is small in the UV and grows larger as we evolve it 
towards the infrared. The sign of this effect is the same as for the non-abelian 
gauge coupling in Yang-Mills theory and it is opposite to the one in QED. The 
main difference is that G k depends quadratically on k while, to lowest order, 
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the gauge coupling in Yang-Mills theory runs only logarithmically. We see that 
gravity is "antiscreening" in the sense that at large distances Newton's constant 
is larger than at small distances. This confirms the intuitive picture that the 
gravitational charge (mass) is not screened by quantum fluctuations but rather 
receives an additional positive contribution from the virtual particles surrounding 
it. 

Let us consider a gravitational (thought) experiment which involves a typical 
length scale r, the distance of two heavy test particles, for instance. If r = k^ 1 
acts as the effective IR cutoff scale, eq.( p.!8| ) suggests the following form of a 
distance-dependent Newton's constant (with factors of h and c restored): 



G{r) = G(oo) 



Gh /l 

0J — + O[-r 



r 2 c 3 



j- j 



(5.21) 



We expectQ that, to leading order in 1/r, the quantum corrected static Newtonian 
potential of two test masses should be closely related to V(r) = —G(r)mim2/r. 
It is interesting to compare ( |5.21 ) to what is actually obtained by a diagrammatic 
calculation of the lowest order correction to the potential. Recently Donoghue ||28f 



has pointed out that quantized Einstein gravity makes a well defined prediction 
for this quantity which is unaffected by the nonrenormalizability of the theory. 
One finds a result of the form 



V(r) 



-G 



17111712 



r 



1 - 



UJ 



Gh 

r 2 c 3 



(5.22) 



The term proportional to (mi + m 2 )/r is a kinematic effect of classical general rel- 
ativity; it is independent of H and is not related to the (3— function of Gk therefore. 
However, the last term in ( |5.22| ), proportional to Gh/r 2 , has precisely the same 
structure as ( [5.21 ). The most recent calculation of uj was performed in ref . [f29[ 
with the result 

-118 , 
u= > 5.23 

157T 

This number has the same sign and is of the same order of magnitude as the value 
found originally in ref. ||28j| , but there is no precise agreement yet. In ref. pj, uj was 



calculated using different methods [[IT], |32| and a negative value was found; this 

5 Recall that in QED the analogous substitution e 2 /r — > e 2 (r~ 1 )/r correctly reproduces the 
leading term of the Uehling potential if the one-loop formula for the running coupling e 2 (p) is 
used l|27tl. 
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would correspond to "screening" rather than "antiscreening" . Possible reasons 
for this discrepancy were discussed in ref.[^9|. While the issue is not fully settled 
yet, it is believed that by correctly identifying and evaluating the set of relevant 
Feynman diagrams, quantum Einstein gravity gives rise to an unambiguous value 
for u. From our investigation of the renormalization group flow we expect this 
value to be positive. 

One can use the full nonperturbative information contained in ( [4.41 ) in order 
to extend the domain of validity of our result towards larger values of gk or smaller 
distances r. This would involve a numerical solution of eq. (|4.38 ) on which we shall 
not embark at this point. 

In our approach we can study the influence of the cosmological constant on 
the running of Gk- It is an interesting question, for instance, whether a large A& 
can destroy the antiscreening character of the gravitational interaction < 0). 
Let us look at ( |5.17| ) with Bi(Xk) given in ( p. 14 ). If there exists a regime with 
i]n > (screening) then -Bi(A) must be positive there. This can only happen if 
the term 5$J(— 2A&) in the brackets on the RHS of ( |5.14j ) is larger then the sum of 
the other terms because the <3>'s are always positive. However, decreases for 

increasing w and finally vanishes for w — > oo. Therefore a negative cosmological 
constant will not change the sign of t]n since B\(Xk) < for A& < 0. 



For Afc > 0, the $'s in (|5.14j) are evaluated at negative arguments w = — 2A&. 
From ( 4.32|) it is clear that blows up for w — > — 1. (The function z + R^(z) 

assumes its minimum value 1 at z = and increases monotonically for z > 0.) 
This signals that our approximation breaks down for A& ~ 1/2 or A& ~ k 2 /2. For 
moderately large values of Afc, Bi(Xk) is still negative. As A& approaches 1/2 from 
below, only the first two terms on the RHS of ( 5.14|) are important. It might be 
that B\ turns negative then, but this would be in a regime where our truncation 
is not reliable any more, and the sign would even depend on R^> in general. 

At this point a general remark concerning the domain of validity of our 
truncation might be in order. In section 3 we showed that truncations of the form 
( 3. 11| ) with Tfc = are consistent with the modified Ward identities provided 
is small. For the Einstein-Hilbert truncation we can evaluate the traces in ( |3.9p 
and we can express Yj, in terms of the functions $^(u>). It is clear, therefore, that 
Yf. becomes large for w — > — 1, and that our truncation cannot account for this 
regime. 
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The running of the (dimensionful) cosmological constant itself is governed 
by the equation 

dAk = VN^k + ^-k 4 G k [l0<^(-2 VA; 2 ) - 8<^(0) - 5r 7iV ^|(-2A A; /A: 2 )l (5.24) 
Ztt 1 J 

If we switch off the renormalization group improvement for a moment and set 

rjN = 0, Afc = on the RHS of eq.( |5.24j ), it has the solution 

X k = -L$i(0)G(A; 4 - A 4 ) + A A (5.25) 

We observe the canonical scale dependence X k ~ k 4 which one expects in any 
naive one-loop calculation: if X k starts off positive at k — A, its absolute value 
decreases when k is lowered until it reaches zero and then X k becomes negative 
(for A large enough). It is obvious that any attempt to fine-tune Aa in such a 
way that lim X k = cannot have a universal meaning because &l{0) depends 
on the form of the cutoff. The evolution equation ( |5.24|) improves on the one- 
loop result in two respects: it includes the effect of the running G k , and via the 
"threshold function" $2 it describes the backreaction of the changing X k on its 
/^-function. In particular, for < and k 2 <C |A^| the relevant IR cutoff in 
the graviton propagator is \X k \ rather than k 2 . Then the graviton modes do not 
contribute to the running of X k any longer, and their decoupling is described by 
the function §l(w). If, on the other hand, the evolution starts with A^ > 0, the 
threshold functions make the coefficient of the /c 4 -term in Q5.24p even larger, and 
the running towards zero is faster than in ( |5.25|) . This effect is counteracted by the 
term r]^X k which is negative for r/jy < 0. It cannot prevent A& from overshooting 
zero, however. 



6 Conclusion 

In this paper we proposed a general framework for the treatment of quantum 
gravity along the lines of the Wilsonian renormalization group. We introduced a 
scale-dependent effective action and we derived an exact renormalization group 
equation which describes its dependence on the built-in infrared cutoff. The ef- 
fective action is invariant under general coordinate transformations; no symmetry 
violating terms are generated during the evolution. It satisfies a set of modified 
gravitational Ward identities which ensure that, in the limit of a vanishing cutoff, 
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the conventional Ward identities are recovered. By virtue of the diffeomorphism- 
invariance of the effective action, fairly simple invariant truncations of the space 
of actions are sufficient to describe the essential physics in a nonperturbative way. 
The modified Ward identities provide a check for the quality of the truncations. 
The evolution equation can be used both for the quantization of fundamental the- 
ories (A — > oo) and for the evolution of effective theories (A finite). It is defined in 
terms of manifestly finite, ultraviolet convergent functional traces. The evolution 
equation by itself is meaningful even if the action is not positive definite. In this 
case the original euclidean functional integral formulation might be problematic, 
and the precise relation between the two approaches is not entirely clear yet. 

As a first application, we have tested our method within a simple truncation 
which retains only the invariants / y/gR and / Jg. Nevertheless, the resulting 
evolution equations for Newton's constant and the cosmological constant contain 
nonperturbative information. In 2 + e dimensions we found corrections to the (3- 
function for Gf. and we determined its dependence on the cosmological constant. 
In 4 dimensions we saw that the /3-function for Gk depends on k quadratically, 
and that Newton's constant increases at large distances. Within its restricted 
domain of validity, this result confirms earlier speculations by Polyakov on a 
possible gravitational antiscreening. 

It would be interesting to allow for a more general truncation and to include 
more complicated invariants in the ansatz for IV Not only higher powers of 
the curvature should be kept but also, and perhaps more importantly, nonlocal 
terms must be included (similar to the 2D induced gravity action / RD~ 2 R, for 
instance). This would lead to a better understanding of quantum gravity in the 
extreme infrared, and might help to clarify certain issues in quantum cosmology. 
For instance, it has been proposed that quantum gravitational effects at large 
distances should be important both in the context of the dark matter problem 
3"1 | and the cosmological constant problem [TJJ, In fact, it is quite clear that 
the nature of the IR divergences, and hence of the renormalization group flow 



for k — > 0, is quite different depending on whether A is zero or not fL3| . In a 
perturbative expansion, one of the traces on the RHS of the evolution equation 
consists of graviton loops attached to external graviton lines. The most singular 
(for k — > 0) diagrams are those which involve the vertices obtained by expanding 
A / y/g, because they do not contain any momentum factors. Hence for A 7^ the 
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renormalization effects should be much stronger than for A = 0, and this could 
eventually drive the cosmological constant to zero. We hope to come back to this 
point elsewhere. 
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